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ABSTRACT 


Vibration control has always been a challenging problem to the helicopter designer. 
The present study addresses the problem on the formulation and solution of an active 
vibration control scheme in helicopters, based on the new concept of Active Control of 
Structural Response (ACSR). First, using a mathematical procedure employing Fisher 
Information Matrix, the optimum sensor locations for vibration measurement in a three 
dimensional flexible fuselage have been identified. It is observed that irrespective of the 
excitation frequency, these optimally selected sensor locations experience relatively high 
levels of vibration. Then, using the measurement from these optimal sensor locations, a 
Multi-Input-Multi-Output (MIMO) control problem has been formulated and solved to ob- 
tain the active control forces required for vibration minimisation in the helicopter fuselage. 
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Chapter 1 

INTRODUCTION 

During forward flight, the time dependent aerodynamic and inertia loads acting on the 
rotor blades cause vibration in helicopters. These rotor loads are transmitted to different 
parts of the fuselage through a complicated load path and cause discomfort to pilot and 
crew equipment deterioration, fatigue damage to the structure and increased maintenance 
cost As a result, these vibratory loads restrict the operation and efficiency 0 f the vehicle. 
With increasing demand for high speed and high performance helicopters, vibration control 
has become an important objective in the design of modern helicopters. References 1-3 
provide excellent review of helicopter vibration and its control. 

Over the years, the vibratory levels in the fuselage of the helicopters have been re- 
duced by using passive vibration control devices and/ or by suitable structural design. For 
present day helicopters, the general requirement is to have a maximum vibratory level of 
0 lg in the fuselage. However, in future, with the adoption of stringent vibration control, 
it will become necessary to reduce the vibratory levels below 0.05g or even 0.02g (Ref.4). 

, a substantial amount of research and development effort has been made 

In recent yeaio, 

at vibration prediction and reduction methodologies in helicopters and also at improving 
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the fundamental understanding of the dynamics of the coupled rotor - fuselage helicopter 
systems. 

The various vibration control schemes adopted in helicopters can be classified as 
either passive or active control methodologies. Since the major cause of vibration in heli- 
copters is the main rotor, one of the passive methods of reducing vibration is to carefully 
design the rotor hub and blades (Ref. 5). Since the rotor design is based on a compromise 
between a number of conflicting aeroelastic, aeromechanical stabilities, performance and 
handling quality requirements, vibration reduction solely by a proper design of rotor blade 
would be impossible. The vibratory levels in helicopters can also be reduced by suitable 
structural modifications (Refs.6, 7) and/ or by employing structural optimisation tech- 
niques (Ref. 8). 

Hub or blade mounted passive pendulum type vibration absorbers (Refs.9, 10) have 
been successful in reducing the vibratory levels in helicopters. But their incorporation al- 
ways leads to increase in rotor weight and aerodynamic drag. The other type of passive sys- 
tem is based on anti-resonant principle. Dynamic Anti-resonant Vibration Isolation(DAVI), 
Anti-Resonant Isolation System (ARIS) and Liquid Inertia Vibration Eliminator (LIVE) fall 
under this category (Refs. 11-14). These isolation devices are usually mounted between the 
fuselage and the rotor or gearbox support system, as shown in Fig.l. They are tuned 
to provide a maximum vibration reduction at a specific frequency. Hence, for any change 
in the operating conditions, there will be a degradation in the performance of these devices. 

Active control methodologies include Higher Harmonic Control (HHC) (Refs.15- 
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19), Individual Blade Control (IBC) (Ref.20), Active Flap Control (AFC) (Ref.21) and 
Active Control of Structural Response (ACSR) (Refs.22-27). It is important to note that 
HHC, IBC and AFC control schemes are provided in the rotating frame, while ACSR is 
employed in the nonrotating frame. Flight tests have shown that higher harmonic pitch 
control of the blade provides a substantial reduction of about 80% — 90% in the vibratory 
levels in the fuselage. However, this control scheme has many drawbacks, such as, (i) rotor 
blade stall at high speeds, (ii) high power requirements and (iii) difficulties in meeting the 
airworthiness criteria. Some of these drawbacks can be eliminated by using an Active Flap 
Control scheme of vibration reduction. Currently, research efforts are aimed at developing 
a smart actuator for Active Flap Control. 

Recently, another vibration control scheme, based on a new concept of Active Con- 
trol of Structural Response (ACSR) has been developed. The idea of ACSR scheme is based 
on the principle of superposition of two independent responses of a linear system such that 
the total response is zero. In the case of helicopters, the fuselage is excited by the appli- 
cation of controlled external actuators at selected locations such that the total response of 
the fuselage due to rotor loads and the external actuator forces is a minimum (Refs.22, 23). 
A schematic of the helicopter system with ACSR scheme is shown in Fig.l. The rotor loads 
(Fhx ,Fhy, Fhzi Mhx, Mry, M H z) are transmitted to the fuselage through the gearbox 
support structure. The support structure is idealised as a spring, damper mechanism and 
a control force generator. In passive scheme, the control force generator corresponds to a 
vibration absorber mass (as in ARIS), whereas in the case of ACSR, the control force gen- 
erator is an active electro-hydraulic force actuator. Preliminary studies based on extensive 
ground and flight tests have shown promising results in reducing vibration in helicopters. 
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The major advantages of ACSR scheme are: (i) less power requirements, (ii) minimal air- 
worthiness requirements because this scheme is independent of the primary flight control 
systems, and (iii) selectively minimise vibratory levels at any set of chosen locations in the 
fuselage. 

The first step in implementing any active vibration control scheme is the mea- 
surement of vibration. In general, the vibratory levels are measured at tail boom/tail 
rotor transmission, cockpit instrument mountings, cabin floor and pilot location (Refe.24- 
26). Even though, these locations may be sensitive, in the light of recent developments 
(Refs.28-30) on the optimal placement of sensors for system identification, an interesting 
question arises: ie., whether the measurement of vibration at the above mentioned locations 
truly represents the vibratory levels in the structure or not. In other words, whether the 
control of vibration at some selected sensitive points in the structure truly corresponds to 
a reduction of vibration in the whole structure or not. A review on the sensor placement 
in distributed parameter(continuous) systems can be found in Ref.31. It is pointed out in 
Ref.30 that the measurement locations play a major role on the quality of measurement 
and in some situation modes may be completely missed. For a simple one dimensional 
structure, the measurement locations can be selected based on experience, but for compli- 
cated three dimensional structures the choice is very difficult. Therefore there is need for 
systematic approach based on mathematical principles to arrive at the optimal sensor loca- 
tions. In Ref.28, Kammer has described a suboptimal procedure for identifying the sensor 
locations in large structures for the measurement of frequencies and mode shapes which 
can be compared with FEM results for correlation studies. This procedure is based on 
using Fisher Information Matrix and Effective Independence Distribution Vector (EIDV) 
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to eliminate sequentially the redundant sensor locations from an initial set of many can- 
didate sensor locations. In Ref.29, this approach was slightly modified, by considering the 
controllability and observability matrices of the system, to identify the actuator/sensor 
placement in a truss structure for modal parameter (natural frequency and mode shape) 
identification. A comparative analysis of EIDV method and Guyan reduction approach is 
presented in Ref.30. The comparative study was based on identifying the sensor locations 
for modal testing of one-dimensional beams and two-dimensional plates. The application 
of EIDV approach for active control of vibrations in helicopters will be highly useful from 
the point of view of practical considerations and such a study has not been reported so far 
in the open literature. 

Coupled rotor - fuselage vibration analysis has been studied by several researchers. 
Though, these studies provide significant contribution towards the understanding of the mech- 
anism of helicopter vibration, most of the studies assume either a rigid fuselage or a 2-D 
beam type fuselage (Refs.14, 32-36). Since the frequencies of the complex helicopter struc- 
ture with heavy concentrated masses will be closely placed (Ref.37), an analysis with 
simplified model of the structure cannot represent all the critical vibration modes of the 
fuselage. Hence, a meaningful study on analysis and control of vibration in helicopter re- 
quires a more realistic structural dynamic model of helicopter fuselage. Therefore, recent 
studies on vibration analysis use a 3-D model for the helicopter fuselage (Refs.27, 38). A 
three dimensional finite element model of a helicopter fuselage, shown in Fig.2, has been 
developed by Mangalick, et. al (Ref.39). The structural dynamic characteristics of this 
model resemble those of a realistic helicopter (Ref.39) and this model is used in the 
present study on vibration analysis and control. 



1.1 Objectives 

The main objectives of the present study are: 


• Identification of optimal sensor locations for measurement of vibration in a 3-D finite 
element model of a helicopter fuselage for active control studies. 

• Formulation of the equations of motion of a coupled gearbox-flexible fuselage heli- 
copter system for vibration analysis. 

• Formulation of an open-loop control scheme for vibration minimisation. 

• Analyse the effectiveness of vibration control using the measurements from optimally 
placed sensors in comparison to the control of vibration using measurements from 
arbitrarily placed sensor locations. 



Chapter 2 

EQUATIONS OF MOTION FOR 
COUPLED 

ROTOR-GEARBOX-FLEXIBLE 
FUSELAGE SYSTEMS 


The analytical study on helicopter vibration requires the development of the equations 
of motion representing the dynamics of the coupled rotor-gearbox-fuselage systems. The 
dynamic model must consist of 

• Rotor blade model 

• Gearbox model 

• Fuselage model 

• Rotor-Gearbox-Fuselage interface model 
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The rotor blade model requires the structural, inertial and aerodynamic operators asso- 
ciated with the rotor blade motion. The flexible helicopter fuselage can be idealised as a 
three dimensional finite element model. The rotor-gearbox-fuselage interface must repre- 
sent the geometry of the interface as well as the aerodynamic interaction in an apporiate 
manner. A schematic of a helicopter system is shown in Fig.l. The main rotor system is 
attached to a rigid shaft from the gearbox. The gearbox is connected to the roof of the 
fuselage through a support structure. 


2.1 Equations of Motion 

A detailed description of the equations of motion of the coupled rotor-gearbox-fuselage 
system is provided below. 

2.1.1 Blade Model 

The rotor blade can be modelled as an elastic beam undergoing deformations in axial-flap- 
lag and torsional modes. The equations of motion of the blade in modal space can be 
written symbolically as 


[AfMy} + + [K] b {y} + {F} = 0 (2.1) 

In this equation, the forcing function {F} contains all the non-linear terms, the aerody- 
namic loading and the effects due to translational and rotational rigid body motions of 
hub. The vector {y} represents the generalised co-ordinates associated with the various 
blade modes. The blade loads contain all harmonics of the rotor rotational frequency. The 
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individual blade loads combine at the hub resulting in the concellation of some of the har- 
monics and the addition of those harmonics which are integer multiples of blade passage 
frequency. The hub loads in the non-rotating frame are obtained after summing up the 
loads associated with the individual blades in the rotating frame. 

Let F k xR ,F k yR ,F k zR and M k xR , M k yR , M k zR represent the hub shears and hub 
moments due to the k — th blade, in the rotating frame, respectively. (The hub-fixed 
non-rotating and the hub-fixed rotating co-ordinate systems are shown in Fig.3). The hub 
loads in the non-rotating frame can be written as 
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(2.3) 


Where NB is the number of blades in the rotor system. The quantities F Hx , F Hy , F Hz and 
M Ri , M Hy , Mjiz represent the hub shears and hub moments respectively. 
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2.1.2 Gearbox Equations 

The gearbox is assumed to be rigid. The gearbox mounting have finite length and each one 
is idealised as a linear spring and a viscous damper element. In addition, a control force 
generating mechanism is acting prarallel to the spring and damper elements. A schematic 
of the mounting is shown in Fig.l. 


The loads acting on the gearbox can be decomposed into two parts, one component 
consisting of the integrated rotor blade loads acting at the hub and the second part repre- 
senting the load from the mountings. While the integrated rotor loads at the hub consist of 
forces and moments, the mountings are assumed to transfer only forces. The force applied 
by each mounting has three components, namely the damping force, the spring force and 
the control force. The gearbox support forces acting at point /, (Fig.l), can be expressed as 


F xi 
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The subscript i indicates the f-th mounting of the gearbox and NC is the total 
number of gearbox supports. I , m and n are direction cosines . of the i-th mounting 
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structure in equilibrium position. The vectors {x u y it Zi} GB and {x 1 ,yuZ i } F represent the 
displacements of the points I and i ' . Points I and /' are the points of attachment of 
the i-th isolator to the gearbox and the fuselage respectively. K x and C { represent the 
stiffness and damping of the i-th gearbox mounting respectively. The control force vector 
{F X i, F y i,F Z i}c can be due to an active isolation device (ACSR) or due to a passive isolation 

device (ARIS). 

The total force on the gearbox, consisting of forces due to isolators, rotor hub loads and 


gravity, can be written as 


\ 


f \ 


f \ 


r \ 

F ' x 
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(2.5) 


The total moment about the centre of mass of the gearbox is given as 
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The rigid body perturbational equations of motion of gearbox are obtained after subtracting 
the mean values of the forces and moments from the total loads. These equations can be 

written as, 

Perturbational Translation Equations 


F x F X m 


m 0 0 


XCB j 

| Fy — Fym 


0 m 0 
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GB 
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(2.7) 
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Perturbational Rotation Equations 

{M — M m }cB = Hqb + & x Hqb 

where, 



C3 = @xGB ® + @yGB j + @zGB k 


( 2 . 8 ) 


2.1.3 Fuselage Equations 

The fuselage is idealised as a three dimensional finite element model, shown in Fig.2. 
The finite elements considered are simple one dimensional beam elements with uniform 
properties. Rotary inertia and shear deformation have been neglected. Lumped masses 
representing equipment, engines, fuel, tail gearbox and end plates are attached to the 
structure at appropriate node locations. The forces acting on the fuselage are those acting 
at the gearbox mountings, tail rotor loads, gravity and fuselage aerodynamic loads. 

The total force on the helicopter fuselage can be expressed as 
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The moment about the centre of gravity of the fuselage can be written as 
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For equilibrium condition, the mean values of the gearbox support forces and mo- 
ments acting on the fuselage must balance the steady components of the fuselage aerody- 
namic, gravity and tail rotor loads. In this formulation, it is assumed that the fuselage 
aerodynamic loads and tail rotor loads are constants. The equations of motion of the fuse- 
lage are as follows 

Perturbational Translation Equations 
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Perturbational Rotation Equations 

{M — ~ HpFOJ'XHp 
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where, 
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(2.13) 


Fuselage Elastic Mode Equations 

The vibratory loads acting on the fuselage are only due to the main rotor. These loads 
are transmitted to the fuselage at the gearbox support points. The fuselage flexible mode 
equations in modal space can be written as, 


[%1 + [<?]{!?} + [JflW = Wf} (2.14) 


where, 

{77} represents the modal vector, 

and Q j F is the generalised force for the j-th vibrational mode, which is given as, 
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Chapter 3 


REDUCED SET OF EQUATIONS 
FOR VIBRATION AND CONTROL 
ANALYS ES 


The coupled set of e*qusai oiwis of rmotuorn fo *itfce ro lor-gearbox-fuselage dynamics, formulated 
in Chapter 2, are simpEibfed for the impose of Ye present study on active vibration control 
in helicopters. In the simplified model, threrotooi ‘blade dynamics are not included. However, 
the hub loads due t- o t ohes e xo tor blades ar*« assuuroed to be acting at the top of the gearbox, 
as shown Fig.4. Se venal assumpYicw lawe Te*ei made in formulating the simplified set of 
equations . 

3.1 AssuumpiL'ions 

1. The gearbox is sassuined to be rignd ^u*dTUndergoes vertical translation (z GB )-, pitch 
{6 y gb) amd roll •(t^ar) rnoiioxs. 
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2. The fuselage is assumed to be undergoing rigid body vertical translation (zp), pitch 
(i 9 y j ?) aid roll {& x f) motions, as well as flexible deformation due to elastic modes. 

3. Only t»tes vertical component, (t), of the hub loads is assumed to be acting at the 
top of the gearbox, simulating a, ground test condition. 

4. The gea-ubox supports are assua-ed to be uni-axial members providing forces, only in 
the z-dmection. 

5. The centre of mass of the gearbox is assumed to be above the centre of mass of 
fuselage on same vertical axis. 

6. The rigiid body rotational motions of the gearbox and fuselage are assumed to be 
small. Hence, the nonlinear te:rms involving products of rotational degrees of freedom 
have hesen neglected. 

7. Th& products of inertia of the gearbox and the fuselage are assumed to be zero. 

3.2 Equations of Motiion of Coupled Gearbox-Fuselage 
System 

The equal io ms of motion of the coupled gearbox-fuselage system can be written in three 
sets. Set 1 describes the rigid body equations of motion of the gearbox; Set II presents 
the rigid body equations of motion of the fuselage and Set III represents the equations of 
motion of the elastic modes of the fuselage. The details of the derivation are given below. 
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3.2.1 Set I: Gearbox Rigid Body Translational and Rotational 
Equations 

The force F { due to tithe spring and damper (Fig.l), acting at the i-th gearbox support can 
be expressed as, 


Ft = Ki {zgb —OyGB Pxi + 6xGB Pyi ~ Zi) + Ci(z G B ~ QyGB Pxi + ®xGB Pyi ~ h) 0(3.1) 

where Z{ and represent the total displacement and total velocity at i-th gearbox support 
point in the fuselage respectively. These can be written as 

Nm 

Zi = zp — 9yp h xi + 0 x p hyi + ^ ( 3 - 2 ) 

3=1 

Nm 

Z{ = zp — 9yp h X { 4 * 8 x p hyi + ^ jz Vj{t) - ( 3 . 3 ) 

3 = 1 

Translational Elation 

NC 

m aB Sgb = F z (t ) - £(F, - F c ‘) (3.4) 

i=l 

Fj represents the control force at the i-th gearbox support. 

Rotational Equations 

nc 

Pitch Equation hvGB ^sGB = Pxi ~ F c ' p xi ) (3.5) 

i=l 

NC 

J-xxGB @xGB — ) ', (Fj Pyi — F c Py%) 

i=l 


Roll! Equation 


(3.6) 
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3.2.2 Set IX: Fuselage Rigid Body Translational and Rotational 
Equations 
'Translational Equation 


wc 

rrifr zp = — Ed) 

i*=l 

Rotational Equations 


Pitch Equation 


Roll Equation 


NC 

Tyyphr = -i(Fih*i-F C ' h *i) 
1=1 


NC 
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(3.7) 


(3.8) 

(3.9) 


3.2.3 Set III: Fuselage Elastic Mode Equations 

The equations of motion of the elastic inodes of the fuselage can be written symbolically 
as, 

[MRij} + [CM +■ [KM = (Or) (3 ' 10) 

The generalised force in the j-th mode, can be written as 


Q i P = j=1 ' 2 (3 ' u) 

Where Nm is the number of elastic modes of the fuselage considered for vibration analysis. 


Equations (3.4-3.10) constitute the complete set of linear, coupled differential equa- 
tions used for the vibration analysis. These equations arenoudimensionalised using radius 
(R) of the blade for length units, mass of blade (m B ) for mass units and rotor speed of 
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rotation (ft) for time. For mathematical convenience, the nondimensional equations are 
written in state space form as 

{ q } = I A ]{?} + [ B ]{[/} + {/} (3.12) 

where, 
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F\ 

F\ 

> 

F i c 

' NCx 1 

The details of system matrices [A] and [B] are given in Appendix A. It may be noted 
that the nondimensional quantities are represented by same symbols as the dimensional 
quantities in Eq.(3.12) as well as in the following sections, for convenience. 




Chapter 4 


SELECTION OF SENSOR 
LOCATIONS FOR VIBRATION 
MEASUREMENT 


The first step in implementing any active vibration control scheme is the measurement 
of vibration. The sensors must be placed at appropriate locations which do not coincide 
with the nodal points of the vibratory modes. For complex three dimensional structure, 
the selection of sensor locations is a difficult, task. In practicej optimal sensor locations for 
vibration measurement are identified by performing extensive shake tests, Refs. (23, 24). 
Therefore, there is strong need for a systematic approach based on mathematical princi- 
ples to arrive at the optimal sensor locations for vibration measurement and control. In 
recent publications, Refs. (28-31), a sub-optimal procedure has been developed for identi- 
fying the sensor locations in large space structures for the measurement of frequencies and 
mode shapes. This procedure is based on using Fisher Information Matrix and Effective 
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Independence Distribution Vector (EIDV) to eliminate sequentially the redundant sensor 
locations from an initial set of many candidate sensor locations. The application of EIDV 
approach for active control of vibration in helicopters will be highly useful from practical 
considerations. 

This Chapter describes the EIDV approach to arrive at the appropriate sensor 
locations for vibration measurement and control. 


4.1 Mathematical Model 

The equations of motion of a flexible structure in finite element domain can be written as 

• [M]f{x} + [C]f{x} + [ K] f {x } = {E} (4.1) 

where [M]p, [ C\p , [K]f and {F} represent the mass, damping, stiffness matrices (of dimen- 
sion NxN) and the external force vector(of size Nxl) respectively. Considering the first 
Nm undamped modes, the modal transformation relation can be written as 

{x} = [$]{? 7 } (4.2) 

Substituting Eq.(4.2) in Eq.(4.1) and premultiplying by [$] T yeilds 

mm + [cm + [Rm = iqf} ( 4 . 3 ) 

where [M], [C], [K] have a dimension of NmxNm and {7?}, { Q F } are vectors of size Nmxl. 

Assuming harmonic input excitation {F} = {F}e wt where a; is a constant, the 
steady state displacement at any point on the structure can be expressed as 
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{?/} = l$s){v} (4-4) 

where the dimension of {y} is M L x 1 and Ml represents the initial number of candidate 
sensor locations. $ s represents the modal matrix corresponding to the initial set of candi- 
date sensor locations. 

To start with, it is assumed that the initial number of candidate sensor locations is 
greater than the number of modal co-ordinates ( i.e.,ML >Nm ). In state feed back control, 
an estimate of the states of the system is required and the best estimate can be obtained 
from the following equations. 

{«} = [S/SJ'VM (4.5) 

V ' 

The underbraced term in Eq.(4.5) is denoted as Moore-Penrose inverse or pseudo-inverse 
(Ref.40) of $ 5 . Since the dimension of 4 > s is (MlxNtti) and Ml > Nm, the rank of 
<3> s is equal to Nm which is same as the number of modal co-ordinates. Hence, there 
are (M L - Nm) rows in $ s which are linearly dependent on the remaining Nm rows. 

Physically, it means that there are {Ml - Nm) additional sensors providing redundant 

information about the Nm modal co-ordinates. From the point of view of reachability of 
certain locations and also due to the cost of sensors, it is not possible to have sensors at 
all locations. Generally, the number of available sensors ( M A ) is less than the number of 
initial candidate measurement locations and it can be greater than or equal to the number 
of modes (i.e., Nm<MA< Ml)- The aim is to eliminate those sensors which provide 
redundant information about the system response. In Eq.(4.5), the symmetric matrix 
[$ s T< i> s ] is denoted as Fisher Information Matrix. Premultiplying Eq.(4.5) by 4> s yields 

*.«} = *.[*/*.]' Vw = { 5 } (4.6) 

> ' 
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If is a non-singular square matrix, then the underbraced term in Eq.(4.6) will be a unit 
matrix. For a general case, let the underbraced term in Eq.(4.6) be denoted by the symbol 
E. 


E = * t [$, T $ a ) (4.7) 

Matrix E is an idempotent matrix i.e., E = E 2 and its eigenvalues are either 0 or 1. In 
addition, the trace of the idempotent matrix E is equal to its rank (Ref.40). Hence, the 
diagonal elements of E represent the fractional contribution to the rank of E and the 
smallest diagonal element (say, Eu ) contributes the least to the rank of E. Since the rank 
of E is equal to the rank of $ s , the i-th row of contributes the least to the rank of 
<f> s . Therfore, the i-th row of can be eliminated without influencing its rank. After 
eliminating the i-th row, the modified 4> s having a reduced size is used to compute the 
new E matrix and the process of elimination is repeated. This procedure is carried out 
sequentially until the number of rows of <E> S is equal to the number of available sensors M A . 
The vector formed by the diagonal elements of E is denoted as the Effective Independence 
Distribution Vector (EIDV). 

Since the inverse of Fisher Information Matrix is required to estimate the modal 
vector, (Eq.4.5), it is important to monitor its Condition Number at every iteration. If 
there is a drastic increase in the condition number, then the elimination process has to be 
terminated. It may be noted that the condition number of a square matrix represents the 
sensitivity of its inverse to very small changes in the elements of the matrix (Ref.41). The 
definition of condition number is provided in Appendix B. 


In every iteration, one can eliminate either one row (one sensor) or a group of 
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rows (group of sensors) whose corresponding diagonal elements ( Eu ) are very small in 
comparison to other diagonal elements (as carried out in Refs. 28, 29). In identifying the 
optimal sensor locations, the advantage of group elimination is that it requires less number 
of iterations as compared to single elimination. However, the disadvantage will be that 
there is a likelihood of increasing the condition number of the Fisher Information Matrix. 
This important conclusion is brought out in Chapter 6, while presenting the results of a 
comparative analysis of single elimination and group elimination. 



Chapter 5 


VIBRATION CONTROL 


During steady forward flight, the helicopter experiences periodic hub loads. The predom- 
inant frequency of the hub loads is NB/ rev, where NB is the number of blades in the 
rotor system. These vibratory loads excite the fuselage structure. The vibratory levels in 
the fuselage axe measured by a set of optimally placed sensors. Using the measurements 
from the sensors, an open-loop (Multi-Input-Multi-Output) control scheme is formulated 
to minimise the vibration in the fuselage. 


5.1 Open-Loop Control Formulation For Vibration Re- 
duction 

The equations of motion of the coupled gearbox-fuselage model have been derived in Chap- 
ter 3. These equations can be written in state space form as (Eq.3.12 ) 

{?} = [4{?} + [£]OT+{/} (5.1) 
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The sensors measure vibratory response of the fuselage. The sensor measurements can be 
written as 


m = [CM ( 5 . 2 ) 

The details of [A], [B], [C], { 9 } and {/} are given in Chapter 3 and in Appendix A. 

For harmonic input {/} = {fje™ 1 , the steady state response can be written as, 

{<?} = [A- iuI)~ l [B]{U} + [A - iuir'if} (5.3) 

The vibratory response measured at preselected sensor locations can be written as 

TO = [cw«} 

= [c\,[a - iuir'mm + mi 

= P1{P} + {6} (5.4) 


where, 

[T] = [C\,[A-iwI\-'[B\ 

{ 6 } = {CUA-^I\-'{f} 


(5.5) 

(5.6) 


Formulating a minimization problem as 

min J = {F} r s {F} s w.r.t {U} 

Substituting for {F} s from Eq.(5.4) yields 

min J = U t T t TU + b T TU + U T T T b + b T b 


Differentiating J with respect to {U} yields 

dJ 


d{U} 


= 2T t TU + 2 T T b = 0 


(5.7) 
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From Eq.(5.7), the best estimate of the control vector minimizing the performance index 
can be written as 

U = -[T T T]~ 1 T T b (5.8) 


Substituting {[/} from Eq.(5.8) and using Eq.(5.2), the controlled vibratory response at 
any location in the coupled gearbox-fuselage system can be obtained. It is important to 
note the control force vector {U} is estimated using the vibratory response at only certain 
preselected locations in the system. For example, these preselected locations could be the 
optimally identified locations or they could represent any set of arbitrarylocations. 



Chapter 6 


RESULTS AND DISCUSSION 


Using the dyn ami c model of the coupled gearbox-flexible fuselage system, several studies 
were performed. The results of these studies are presented in three sections. First section 
describes the results pertaining to the study on the choice of sensor locations for vibration 
measurement. The second section presents the results of a vibration analysis of the coupled 
gearbox-fuselage model. The results of the study on vibration control in fuselage are 
presented in the third section. 


6.1 Choice of Sensor Locations for Vibration Mea- 
surement 

Using the mathematical formulation presented in Chapter 4 on EIDV approach, optimal 
sensor locations have been identified in a three dimensional finite element model of a 
helicopter fuselage. Before presenting the results corresponding to the helicopter fuselage 
problem, an example problem has been solved to highlight the essential features of the 
EIDV approach. 
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6.1.1 Example Problem 

Consider an example in which there are six candidate measurement locations and {y u i = 
1 , • • • , 6 ) and three modes represented by the generalised co-ordinates (rji, 772, rja), which is 

given as 

{y} = l *.]{*?} {6A) 

where, 



1 0 0 
0 1 0 
0 0 1 
2 0 0 

3 0 0 

4 2 0 
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The three column vectors of <!>,, correspond to the modal displacements at the six mea- 
surement locations. In this example, the rank of is 3. The 4-th and 5-th rows of 
are integer multiples of 1-st row. The 6-th row is a linear combination of rows 1 and 2. A 
sensor placed at location 6 will measure a maximum displacement in mode 1 and mode 2. 
A sensor at location 3 is essential to measure the 3-rd mode. In this case, a minimum of 
three sensors are required to estimate the three generalised co-ordinates. It is essential to 
have one sensor at location 3, preferable to have the second sensor at locations 6 and the 
third sensor can be placed at anyone of the remaining four (1,2,4 and 5) locations. Using 
the EIDV approach, the redundant sensor locations are eliminated sequentially one at a 
time. The- results of the sequential elimination are provided in Table 1. 

Table 1 contains [$ s ]j and [E] - (subscript j refers to j-th iteration) at every iteration 
along with the lowest diagonal value of [E]-. In the first iteration, the lowest diagonal 
value of E is E\\ which is Therefore the first row has the least influence on the rank 
of Eliminating the first row (i.e., the sensor at location 1), the new [$ s ]2 is formed 
for the second iteration. In the second iteration, the lowest diagonal value of [E] 2 is E& 
which is |f. Eliminating this row (note that this row corresponds to the sensor location 4 
in the original set of six sensor locations), the new [<& s ]3 for the third iteration is formed. 
In the third iteration, the lowest diagonal value of [i?] 3 is Eu which is |f. Eliminating 
this row (i.e., the sensor at location 2 in the original set of six locations), the final [$^4 
is formed. In this case,- [E)\ is an identity matrix, indicating that all diagonal elements 
(i.e., all sensors) are equally important. Further elimination of any one row will render 
the Fisher Information Matrix a singular matrix. This means that we are attempting 
to estimate three generalised (independent) co-ordinates with two sensor measurements, 
which is not possibc. In this example, the condition number of the Fisher Information 
Matrix exhibits a very small variation with iteration. 

The final set of optimal locations for the sensors are at locations 3, 5, and 6. It is 
interesting to note that in all iterations, the value of the diagonal element corresponding 
to the sensor at location 3 is always one, indicating that this sensor is very important as 
noted earlier. The trace of the matrix E in every iteration is equal to 3, which is the rank 
of <I> 5 . Another important observation is that the diagonal elements of E cannot be greater 
than 1. A formal mathematical proof for this statement can be found in Ref.28. 
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6-1.2 Helicopter Fuselage Problem 

Figure 2 shows the baseline finite element model of a helicopter fuselage. The length 
of the helicopter model is 8.25m, the height is 2m, and the width is 3m. The fuselage 
is 4m long, having a width of 2.5m and a height of 1.5m. The tail boom is 4.25m in 
length, with a horizontal stabilizer having a span of 3m attached near the end. The total 
structural mass is 418.6kg. In addition, lumped masses representing two engines (each 
weighing 180kg), tail gearbox (weighing 40kg) and two end plates (each weighing 2kg) are 
also attached to structure at appropriate nodes. For the basleine configuration, the weight 
of one engine is uniformly distributed at nodes 41, 45 and 46; the other engine weight is 
unformly distributed at nodes 44, 48 and 49. The tail gearbox is placed at node 64. The 
weight of end plates are lumped at nodes 56 and 62 of the horizontal stabilizer. Total 
number of nodes and the degrees of freedom of the finite element model are 64 and 384 
respectively. The details of the structural properties, node locations and other data are 
given in Ref.39. It was shown in Ref.39 that the undamped natural frequencies and mode 
shapes of this model are similar to those of a realistic helicopter. This model is used in the 
present study on the selection of sensor locations and vibration control. Figure 5 shows 
the first four modes of the fuselage. Table 2 presents the natural frequencies of the first 20 
modes of the baseline configuration of the helicopter along with the nondimensional values. 

Assuming that the main rotor system consists of four blades, the vibratory hub 
loads will have a nondimensional excitation frequency of 4/rev. It is evident from Table 2 
that, the frequency of the rotor excitation loads is well below the natural frequency of the 
20-th mode which is 6.4064. Therefore, it is more than sufficient to consider the first 20 
modes for the vibration analysis of the helicopter fuselage. (The validity of this statement 
is proved by performing a convergence study on the vibration analysis, the results of which 
will be presented in Section 6.2). 

Considering the first 20 modes of the fuselage (A r m=20), the modal matrix <J> 5 is 
formulated. Since the vibratory level in the vertical (z) direction is more predominant, 
without loss of generality, it is assumed that the sensors measure only the z-component of 
the fuselage vibration. Therefore, in the formulation of the modal displacement in the 
z-direction only is considered. Initially it is assumed that all the 64 nodes are the candidate 
sensor locations (i.e., M L = 64). The redundant sensor locations were eliminated by single 
elimination (one at each iteration) or by group elimination. 
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The results corresponding to single elimination procedure are shown in Figs.6 and 7. 
The variation of condition number of the Fisher Information Matrix with iteration number 
is shown in Fig.6. In this figure, the quantities within the bracket represent the number 
of sensors and the arrows indicate the value of condition number at the corresponding 
iteration. In the first iteration, there are 64 sensors and the condition number for this 
case is 2.86E+05. It is observed that the condition number remains more or less same till 
39-th iteration. Then in the next 6 iterations, the condition number increases from around 
3.55E+05 to 7.22E+05, but the order remains the same. At the end of 44-th iteration, 
there are 20 sensors and the corresponding condition number is 7.22E+05. These optimally 
selected 20 sensor locations are shown by node numbers in Fig. 7. 

The result of group elimination is shown in Fig.8. In this case, in every iteration, 
group of sensors, whose corresponding Ea values in Effective Independent Distribution 
Vector (EIDV) are almost equal, are eliminated. The variation of condition number with 
iteration is shown in Fig.8. In the first iteration, there are 64 sensors and the corresponding 
condition number is 2.86E+05. Till 13-th iteration, the condition number remains more or 
less in the same order. There are 25 sensors in the 13-th iteration and the corresponding 
condition number is 9.02E+05. With further elimination of sensors, the condition number 
increases drastically. In the 14-th iteration, there are 23 sensors with a condition num- 
ber equal to 1.69E+13. Since there is a drastic increase in the condition number, the 
elimination process is terminated at the 14-th iteration. 

Comparing the results shown in Figs.6 and 8, it is interesting to note that even 
with 23 sensors, identified by the group elimination process, the condition number is of the 
order of 1.0E+ 13, whereas with 20 sensors, identified by the single elimination process, the 
condition number is of the order of 1.0E+05. 

6.1.3 Failure Analysis of Sensors 

A preliminary study is carried out to analyse the effect of sensor failure on the quality 
of state estimation. The parameter used for comparing the quality of estimation is the 
condition number of the Fisher Information matrix. It is assumed that initially a redundant 
set of 25 sensors are available for the estimation of response of 20 modes. Two sets of 25 
sensors were chosen; one set is obtained by group elimination and the other set by single 
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elimination. 

The study of failure analysis is carried out for two cases. In the first case, any one 
sensor is assumed to fail at a given instant. It may be noted that there 25 combinations 
of failure for 25 sensors. In the second case, it is assumed that any two sensors can fail 
at a given instant; for this case there are 300 combinations of failure for 25 sensors. For 
all these combinations of failure, the condition number of Fisher Information matrix is 
calculated after eliminating the row, corresponding to the failed sensor, from the matrix 
Then, the condition numbers corresponding to the various combinations of failure are 
arranged in ascending order. Figure 9 shows the variation of condition number for one 
sensor failure analysis and Fig. 10 shows the variation for two sensors failure analysis. The 
arrows in the figures indicate the range of variation of condition number. It can be seen 
from Fig.9 that the range of variation ofcondition number, for sensors obtained by single 
elimination, is very small compared to that for the sensor set obtained by group elimination. 
A similar observation can also be seen in Fig. 10, for two sensors failure analysis. In this 
case(Fig.lO), the range of variation of condition number for group elimination case is of the 
order of 1.0E+08, while the range of variation for the single elimination case, is of the order 
1.0E+01. These results clearly indicate that single elimination approach of selecting the 
optimal sensor locations provides a very stable condition number for Fisher Informatiom 
matrix than the group elimination process. 

6.1.4 Condition Number and Normalisation of Modal matrix 

It can be seen from the results that even for single elimination, the absolute value of the 
condition number (order of 1.0E+05) of Fisher Information matrix is quite high. However, 
in the selection of sensor locations, the variation in the order of condition number is more 
important than the absolute value itself. The reason is that the condtion number of the 
modal matrix can be arbitrarily varied by employing different normalisation factors for 
different modal vectors. This point is illustrated in the following example. 

Consider a matrix $ s , whose size is 3x2. The two modes correspond to the first 
two modes of uniform cantilever beam. A modified matrix <3> s is obtained by multiplying 
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two columns of $ s by different multiplication factors. 



0 0 
0.2776 -0.5287 
1.00 1.00 



0 0 
0.02776 -0.005287 
0.1 0.01 


The condition number of the Fisher Information matrix [$ s r $ s ] is 7, whereas the condition 
number of [<f s $ s ] is 207. It is interesting to note that the condition number of modal 
matrix can be arbitrarily changed by suitable multiplication factors. 

6.1.5 Sensitivity of Sensor Locations to Modifications in Fuselage 
Structure 

A study was conducted to analyse the sensitivity of sensor locations to modifications in 
the finite element model of helicopter fuselage structure. From the baseline configuration, 
the fuselage structure was modified by shifting the weight of one engine from node 48 to 
node 49. All the other data remain same. The natural frequencies and mode shapes of 
the modified structure were calculated. Table 3 presents the natural frequencies of the 
modified structure along with those of the baseline configuration. From this Table, it can 
be seen that there is very little change in the natural frequencies of the structure due to the 
modification. Using the modified structure, optimal locations for 20 sensors were obtained 
by single elimination approach and these locations are shown in Fig.ll. Comparing Figs.7 
and 11, it is noted that from the baseline case, the sensor at node 12 is moved to node 8; 
sensor at node 20 is shifted to node 27 and sensor at node 43 is shifted to node 39. 
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Ihis study indicates that tiie E1DV approach for the selection of optimal sensor 
locations seems to be rather reasonably sensitive to structural modifications. 

6.1.6 Validation of Optimal Sensor Locations 

A vibration analysis was performed using the coupled gearbox-fuselage equations, by ap- 
plying a vibratory force at the top of the gearbox. Total number of degrees of freedom 
considered in this analysis are 26. These include 3 rigid body modes of the gearbox (heave, 
pitch and roll), 3 rigid body modes and 20 flexible modes of the fuselage. The data used 
for the vibrational analysis are given in Table 4. 

The vibratory levels in the fuselage were calculated for different excitation frequency, 
namely, 1/rev, 2/rev, 3/rev, 4/rev and 5/rev. The results of this analysis are shown in 
Figs. 12-16. In these figures, the vibratory levels (g-levels) at different nodes are indicated 
by impulses. The arrows indicate the optimal locations for the sensors obtained by single 
elimination process. It is interesting to note from these figures that irrespective of the 
input excitation frequency, the optimally placed sensors measure high levels of vibration. 
For 1/rev excitation, the sensor at node 64 measures the highest level of vibration of 0.48g 
(Fig.12). For 2/rev excitation, though there is no sensor at node 38, where the g-level is 
maximum, there is a sensor at node 37 measuring the second highest level of vibration of 
0.65g (Fig.13). For 3/rev excitation, the peak response is at node 38. But there is a sensor 
at node 37 measuring the second highest level of vibratory response (Fig. 14). For 4/rev 
(Fig.15) and 5/rev (Fig.16) excitations, the peak response occurs at node 33. But there 
are two sensors at node locations 32 and 34 measuring the second highest level of response. 
These results indicate that the optimally selected sensors measure high levels of vibration. 


6.2 Vibration Analysis 

Assuming that the rotor system consists of four blades, a 4/rev hub vibratory load is 
applied at the top of the gearbox. Since the vertical load is predominant, for the purpose 
of analysis, only vertical component of load ( F z {t) ) is applied. The. gearbox is assumed 
to supported on the roof of the fuselage at four nodes. Set I consists of the four nodes 39, 
48, 46 and 37 and Set II consists of nodes 43, 47, 42 and 38. The data for this model is 
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given in Table 4. 

First, a convergence study on the vibratory levels in the fuselage with respect to the 
number of fuselage modes, was performed. The gearbox is assumed to be attached to the 
fuselage at nodes given in Set I. Figure 17 presents the results of the convergence study. 
Inclusion of 20 flexible modes of the fuselage provides good convergence in the vibratory 
levels at all nodes. The reason for the good convergence is that all the important modes 
whose natural frequencies are in the neighbourhood of the excitation frequency 4/rev are 
included in the analysis. This is evident from Table 5 presenting the natural frequncies 
of the coupled gearbox-fuselage system for the several cases considered in the convergence 
study. Hence, in all subsequent studies on vibration measurement and control, 20 flexible 
modes of fuselage are included. 

The influence of locations of attachment of the gearbox to the fuselage is also 
studied. Set I and Set II represent the two cases of connectivity of the gearbox to the 
fuselage. The vibratory levels at various nodes in the fuselage for these two cases are 
shown in Fig.18. Set I connectivity provides a higher level of vibration in the fuselage than 
Set II. Hence, Set I connectivity is chosen as the reference case for vibration minimisation 
studies described in the following section. 


6.3 Open-Loop Vibration Control 

For the vibration control studies, the total number of degrees of freedom of the dynamic 
model is 26. This consists of 3 rigid body degrees of freedom of the gearbox, 3 rigid body 
degrees of freedom and 20 flexible modes of the fuselage. Since, there are 23 degrees of 
freedom for the fuselage, 23 sensor locations were identified by EIDV approach employing 
single elimination process. These optimally selected 23 sensors locations are indicated by 
node numbers in Fig. 19. 

Incorporating an open-loop control scheme, described in Chapter 5, an attempt is 
made to minimise the vibratory response of the fuselage. The control forces, required for 
vibration minimisation, are evaluated using measurements from several sets of sensor loca- 
tions. These different sets of sensor locations correspond to (i) optimally placed 23 sensors, 
(ii) arbitrarily placed 5 sensors(node locations 12, 13, 20, 21 and 64), (iii) arbitrarily placed 
10 sensors (node locations 1, 2, 20, 21, 31, 35, 50, 56, 62 and 64) and (iv) arbitrarily placed 
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23 sensors (node locations 1 through 22 and 64). The gearbox is attached to the fuse- 
age at node locations given in Set I. Assuming a. 4-bladed rotor system,. the nondimensional 
frequency of the excitation force is 4/rev. The relevant data are given in Table 4. 

Using the g-levels measured at the optimally selected 23 locations, the control forces 
required for minimisation of vibration in fuselage were calculated. Figure 20 shows a 
comparison of the baseline vibratory levels along with the controlled response. The fuselage 
vibratory level has been reduced substantially from the basline peak acceleration of 0.284g 
to a level of 0.5E-04g at node location 33. In addition, the vibratory levels at all nodes in the 
fuselage are reduced to very low levels. But the gearbox C.G. experiences an increase in the 
g-level, i.e., the gearbox g-level increased from a value of 0.0477g to 0.0625g. Figures 21(a) 
and 21(b) show the magnitudes and the phase angles of the four control forces required 
for vibration minimisation. The control forces are almost 180 degrees of out-of-phase to 
the applied external force. This observation can be mathematically proved using a simple 
two degree of freedom model with a control force generator placed between them, which is 
described in Appendix C. 

To analyse the effectiveness of vibration reduction using optimally placed 23 sensors, 
a vibration reduction analysis using different sets of sensors located at arbitrary nodes in 
.the fuselage, was performed. The controlled response for these cases of arbitrary sensor 
locations are compared with the controlled response for the case of 23 optimally placed 
sensors. These results are shown in Figs. 22-24. 

Figure 22 shows the results of the controlled vibratory response, obtained using 5 
arbitrary sensors and 23 optimally placed sensors. It is observed that the peak acceleration 
of the controlled response with 5 arbitrary sensors is 0.31E-02g at node location 43. The 
peak acceleration of the controlled response for optimally placed 23 sensors is 0.138E-03g 
at node location 5. Figure 23 shows a comparison of the controlled vibration with 10 
arbitrary sensors and 23 optimally placed sensors. For the case of 10 arbitrary sensors, 
the peak acceleration of the controlled response is 0.211E-03g at node location 33, which 
is about 53% more than that for the case of optimally placed 23 sensors. Figure 24 
shows the controlled response for the case of 23 arbitrary sensors along with the controlled 
response with 23 optimally placed sensors. For the case of 23 arbitrary sensors, the peak 
acceleration is found to be 0.203E-03g at node location 41, which is 47% more than that for 
the case of optimally placed 23 sensors. The magnitudes and phase angles of the control 
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forces for all these case are shown in Figs. (25-27) and are also presented in Table 6. It is 
interesting to note that though there is very small variation in the magnitudes and phase 
angles of the control forces, there seems to be a large variation in the peak acceleration 
of the controlled vibratory response of the fuselage. These results clearly indicate that 
the vibration control using measurements from the optimally placed sensors provide the 
minimium peak acceleration in the fuselage. This shows that the EIDV approach with 
single elimination seems to be an effective technique for identifying the sensor locations for 
vibration control in helicopters. 

It is observed that in all these vibration minimisation studies, even though there 
is vibration reduction in the fuselage, the gearbox c. g. experiences an increase in the 
acceleration level (Figs. 22- 24). So an attempt was made to reduce simultaneously the vi- 
bratory levels at the gearbox c. g. as well as at the fuselage. In this case, the control 
forces required for vibration minimisation were obtained using measurements from 24 sen- 
sors' (23 optimal sensor locations in the fuselage + 1 sensor at gearbox c. g.). Figure 28 
shows the comparison of controlled vibratory levels obtained by using measurements from 
23 optimal locations and those for the case of 24 sensors. It is interesting to observe that 
with 24 sensors there is no improvement in vibratory levels at the gearbox c. g., but there 
is deterioration in the fuselage vibratory levels. The magnitudes and phase angles of the 
control forces for the 24 sensors case arc shown in Fig.28, and are also given in Table G. 



Chapter 7 


CONCLUDING REMARKS 


The problem of vibration reduction in helicopter fuselage, using a new concept of Active 
Control of Structural Response (ACSR), has been formulated. The equations of motion 
representing the dynamics of a coupled gearbox-fuselage model have been derived. Using 
these equations, several studies have been performed. They are (i) identification of optimal 
sensor locations for vibration measurement, (ii) influence of number of flexible modes on 
the vibratory response of a three dimensional finite element model of a fuselage structure 
and (iii) formulation and solution of a Multi-Input-Multi-Output (MIMO) control scheme 
for vibration minimisation in the helicopter fuselage. 

The important conclusions of this study are summarised below. 

1. A detailed description of the Effective Independence Distribution Vector (EIDV) 
approach for the identification of sensor locations for vibration measurement is pre- 
sented. An example problem has been formulated and solved to highlight the essential 
features of this EIDV approach. 
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2. With condition number of Fisher Information matrix as a reference parameter, it was 
observed that the single elimination process provides a condition number which is 
many orders lower than that corresponding to group elimination process. 

3. The condition number of modal matrix is shown to be highly dependent on the 
normalising factors used for different modal vectors. 

4. Irrespective of the input excitation frequency, the optimally identified sensor locations 
by the single elimination process, experience high levels of vibration. 

5. For an excitation frequency of 4/rev, it was observed that 20 flexible modes of fuselage 
provide a good convergence for the vibratory response of fuselage. 

6. Vibration control using measurements from the optimally selected sensor lcoations 
provide maximum reduction in the g-levels of the fuselage vibration as compared to 
the controlled response using measurements from arbitrarily placed sensor locations. 

7. When the vibratory levels in the fuselage are minimised, the gearbox experiences 
a higher level of vibration in comparison to the baseline g-level. While trying to 
minimise simultaneously the vibratory levels in the fuselage and the gearbox, it was 
observed that there is no reduction in the vibratory levels at gearbox; but there is a 
deterioration in the control of vibration in the fuselage. 



7.1 Scope for Further Work 

This study can be further extended to incorporate 


• Closed-loop active control methodologies for vibration minimisation 

• Inclusion of rotor dynamics in the analysis 


Vibration analysis and minimisation at several flight conditions. 
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Table 1. Example Problem 
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Mode 

Natural frequencies of helicopter model 

Frequencies in Hz 

Nondimensional frequencies 

1 

3.5109 

0.6894 

2 

4.1491 

0.8146 

3 

9.9545 

1.9545 

4 

12.0511 

2.3662 

5 

13.0459 

2.5615 

6 

13.2619 

2.6039 

7 

15.4780 

3.0391 

8 

16.7103 

3.2811 

9 

18.4370 

3.6201 

10 

19.7831 

3.8844 

11 

22.5630 

4.4302 

12 

23.8393 

4.6808 

13 

24.4473 

4.8002 

14 

24.7372 

4.8571 

15 

25.0777 

4.9239 

16 

26.6665 

5.2359 

17 

30.3410 

5.9574 

18 

30.4265 

5.9742 

19 

31.3445 

6.1544 

20 

32.6280 

6.4064 


Table 2: Natural frequencies of first 20 modes of the baseline configuration of helicopter 
fuselage 




Mode 

Natural frequencies 

Baseline fuselage 

(Hz) 

Modified fuselage 

(Hz) 

1 

3.5109 

3.5371 

2 

4.1491 

4.1294 

3 

9.9545 

9.8118 

4 

12.0511 

11.7661 

5 

13.0459 

12.9444 

6 

13.2619 

13.2014 

7 

15.4780 

15.2866 

8 

16.7103 

16.6090 

‘9 

18.4370 

18.2183 

10 

19.7831 

19.9171 

11 

22.5630 

22.5298 

12 

23.8393 

23.9533 

13 

24.4473 

24.4478 

14 

24.7372 

24.5707 

15 

25.0777 

25.0781 

16 

26.6665 

28.1926 

17 

30.3410 

30.3799 

18 

30.4265 

31.1950 

19 

31.3442 

32.3254 

20 

32.6280 

33.9449 


Table ' 3: Natural frequencies of modified helicopter fuselage system along with baseline 


configuration 
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Table .4: Data for vibration analysis and control 
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Mode 

Nondimensional natural frequencies 

No. 

Rigid Body Modes 

RBM + 

RBM + 

RBM + 

RBM + 

RBM + 


( RBM ) 

5 modes 

10 modes 

15 modes 

17 modes 

20 modes 

1 

7.6932 

0.6641 

0.6639 

0.6639 

0.6637 

0.6637 



0.8146 

0.8145 

0.8145 

0.8146 

0.8146 

b 

12.3890 

2.2913 

2.2299 

2.2057 

2.2004 

2.1963 



2.5082 

2.3803 

2.3693 

2.3121 

2.2457 

B 


2.7589 

2.4414 

2.3957 

2.3851 

2.3838 

n 


7.9603 

2.6042 

2.6025 

2.6020 

2.6020 

n 


10.1363 

2.7016 

2.6932 

2.6529 

2.6312 

8 


12.6932 

3.0658 

3.0605 

3.0594 

3.0590 

9 



3.5981 

3.5914 

3.4960 

3.4744 

10 



4.6597 

3.7922 

3.7104 

3.7010 

n 



9.6157 

4.4329 

4.4303 

4.3956 

B 



10.4418 

4.7436 

4.4594 

4.4291 

B 


| 

14.5201 

4.7761 

4.7694 

4.7667 

B 




4.8916 

4.7981 

4.7981 

15 




6.1996 

4.9236 

4.9235 

16 




9.7430 

5.4685 

5.4413 

17 




11.0455 

6.5694 

5.5312 

18 




14.7380 

11.7584 

5.8245 

19 

I 




11.9704 

6.1304 

20 





15.0572 

7.0591 

21 






11.9308 

22 






12.3048 

23 






17.2981 


Table 5: Natural frequencies of the coupled gearbox-fuselage system for several idealisa- 


tions 
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Figure ■ 4: Coupled gearbox-fuselage dynamic model 
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four mode shapes of the fuselage 
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ITERATION NUMBER 

Figure 6: Condition number of Fisher Information Matrix vs iteration num- 
ber. (Single Elimination) 



Figure 7: Optimal sensor locations for baseline configuration of the helicopter 
(20 sensor locations indicated by node numbers) 
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Figure - 8: Condition number of Fisher Information Matrix vs iteration number. (Group 
Elimination) 





NUMBER OF CASES 

Figure - 9: Range of variation of condition number with one sensor failure 
(Initial set of sensors: 25) 



Figure 10: Range of variation of condition number with two sensors failure 
(Initial set of sensors: 25) 






g-level 



Figure 12; Vibratory response of the fuselage for 1/rev hub force excitation 
(Optimal sensor locations indicated by arrows) 





Node Number 

Figure 13: Vibratory response of the fuselage for 2/rev hub force excitation 
(Optimal sensor locations indicated by arrows) 



Figure 14: Vibratory response of the fuselage for 3/rev hub force excitation 
(Optimal sensor locations indicated by arrows) 






Figure • 15: Vibratory response of the fuselage for 4/rev hub force excitation 
(Optimal sensor locations indicated by arrows) 



t fnQPlfltre for 5/rev hub force excitation 
Figure 16: Vibratory response of the fuselage t / 

(Optimal sensor locations indicated by arrows) 
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Node Number 


Figure 18: Comparison of vibratory response for two sets of connectivity of gearbox to 
fuselage (Set I - nodes 39,48,46,37, Set II - nodes 43,47,42,38) 
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Figure ' 20: Comparison of baseline and controlled vibratory levels (Optimal 
sensor locations indicated by arrows) 
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Figure 22: Effectiveness of vibration control with 23 optimal sensors w.r.t 
vibration control with arbitrary 5 sensors 
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Figure 23: Effectiveness of vibration control with 23 optimal sensors w.r.t 
vibration control with arbitrary 10 sensors 



Node Number 


Gearbox 

^C.G. 


"i i i i i i i i i i i — i — i — i — \ — i — i — i — r 

Vibration Control with optimal 23 sensors* 
Vibration Control with arbitrary 5 sensors □ 


IWp 


□ □ 


w. 




? tp 

toT 


% 


W 




TO 




W 00 


m 


■it tff 


1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 
Node Number 





73 


o.i 
0.01 
0.001 

1 

i 0.0001 

le-05 

le-06 

le-07 

Figure 24: Effectiveness of vibration control with 23 optimal sensors w.r.t vibration 
control with arbitrary 23 sensors 
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Figure 25: Magnitude and phase angle of control forces for abritrary 5 sensors 
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Figure 26: Magnitude and phase angle of control forces for abritrary 10 


sensors 








Phase angle 



Node Locations for control forces 

(b) Phase angle 

Figure 27: Magnitude and phase angle of control forces for abritrary 23 


sensors 
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Figure 28: Effectiveness of vibration control with 23 optimal sensors w.r.t vibration 
:ontrol with 24 sensors 
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Appendix A 


System matrices [A],[B] and [C] 


A 


{0} : 

i 

“i : 

0*2 

“3 i 

a4 

as 

a6 


2(6+Nm)x2(6+Nm) 


_ 


* 

a x 

= 

On “12 “13 

0>2 

= 

021 022 “23 

a3 

= 

0 3 i “32 “33 ] 

a4 

= 

041 “42 “43 ] 

a 5 

= 

£ 051 “52 “53 ] 

a6 

= 

| 061 “62 “63 ] 
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{-T^KiPy, 


ES? ifiP» 
-ESS^Si 

iPxiPyi 


-ES jKtfn) 

-eSS^pV} 


mcB 


i=l 
^ NC 


l yvGB 


{-ES«to 


{ES 


IxxGB 


-eSS^ eSS^M 
eS S* iM* -e£?*>M*} 
-EJSSM* ESS ^VP*) 


=*-{Efei«*uXi 

mcs L 1 1 

rM- £*5 »*>**<»* 

*WGB 


Ei=l^ 2 zKi ••• 

-E^^z^Pxi ••• - T.iJi *Vmz^P«} 
^2z KiPyi • • • EiSl ^ Vmz KiPyi } 


3xNm 
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i r v-JVC n. 

— i 2^=1 ^ 

1 f x^NC , 


1 VVGB 


{E£, Op* 


{-ESC,?; 


IxxGB 


y i 


El C iCi Pxi 

-H^CiV 2 xi 

Y,?= C lCiP X iPyi 


-El^CiPyi} 
E ^CiPyiPti} 

-e£W*} 





m G B 


{” SiS CiPxi 


1 yy G B 


{ES CiPy i 


IxxGB 


— Y^i=i Cih x i Eila Cihyi) 

Efe* OihxiPxi — Efcl C{hyiP X i\ 

— Eil'l CihxiPyi Ei^l CihyiPyi} 


a 23 





mcB 

L - {-Effi^l ztfiP* 


l yyGB 

1 


IxxGB 


{ZZZ&uCiPyi 


-E^^zQPxi 
ZiL C l VlzCiPyi ' ' • 


Efc* tfVmzCi} 

- Ei^l ^VmzC'tPxi} 

Ei^ ^’.VnizCiPyt} 


3xNm 
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771 p 


1 wf 


{EK Ki 
{-ES5*A 




-T.liK<p« EfiS^} 

Z&KtoK, -EjSWte} 


— Ei="l KiPxihyi H'i=i KiPyihyi} 


3x3 




Eili Kihzi — E^i Kith*} 

-•ZZSKilf* Egg ft v^} 

ZliKMhyi -ZHgKi AV} 


^{ESS *..**- 


-T.iA^Kihyi 


-ESt^WQ 
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i r 

Ci 

-Z?= C xCip xi 

eS c iP!rf } 


vuf 

C{h X i 

Et c i C iPxi h xi 

-EgOvM 


Tir (EffiCVi* 

-Zl C iC lPxi h yi 

ESS CiPyikyi) 

3x3 




L,i=l 


{EgCA 


l VVF 


{- ES CiA,, 


IxxF 


TS&CtKi -ESfiW 

— Sil? Cih 2 x i Sil? Ci/lyi/irci} 

YSZCM* -El C iCih 2 yi } 



'ki-TZiVuCi 


-ES^i 

ZZZ&lzCihzi 

-E^^Cihyi 


- T*£ ^NmzCi} 
Ysi^l NmzCihxi} 
— f}rnzCihyi\ 


3xNm 
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051 


= m 


-1 


ESl WuKi 

TSS. 


Egg VuKiPxi EiS VlzKiPyi 

■ Effi &2xKiPxi ^l^zKiPyi 


E£l ^Vm/i - EfcS ^NmzKiPxi E ^Nmz^iPyi 


NmyS 


052 


= W 


- EfJi ^uKi Ex^ VuKihxi -TgPiVuKihyi 

- Efcl Effi - Effi y i 2 zKihyi 

-Y^l^NmzKi l£=l VNmzKihxi ~ Ei=l ^Vm Z #A* J ^3 





Effi ^2*^1* 


EiS Ki&u&Nnu Ei=l Ki$f l 2 z^ t Nmz 


'ZgSKi&Nmz&U 
ZSiKi&SmM* 2* 




Nmz 


NmxHm 


-iMr'ifl 
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in-i 


[M] 


ElS 

eK 


■ sS'e'uCiPri sS^LCiJv 

■ EStfj . Ca * E 


Eg*WC, -E Effi^CflV 


NmxZ 



[ M ]- 1 


-EfcS^ixC'i ZgZ&uCih* 

-EfJi^zCi Hi=i^zCih xi 


-JM&uCih# 

-ZfJiVizCihyi 


Ete? ^Vm Z Ci Efel ^NmzCihxi - E,=? ^VmzC^J/i J ^ mx3 
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EfiCifSz ^ 


EiS 


Efe * C ' i ^ lz^Vmz 




Elrf^ATnu^. 

E^C'i^Nmz^Vmz 
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-[ mi - HCi 



86 





i 
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IxxGB 


1 

mcB 

fyVGB 

PyNC 

IxxGB 


3xNC 






^ Nmz 


* 2 lz ■■■ 

^2z ... 

^ 2 Nmz • ‘ - 


* NC u 

* NC 2 z 
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i,i 



: 


(n+l)xl 



.1 

0 

0 

0 

0 

0 

0 

0 

0 : 



0 

0 

0 

1 

—Xpi 

~Vf i 


^2 ••• 

Nm : 

' 

c = 

0 

0 

0 

1 

-Xf 2 

~VF2 

tf 2 ! 

# 2 2 ... 

^ 2 Nm ■ ■ {0} 



0 

0 

0 

1 

~X F n 

—UFn 


\fl Nm 2 . . . 

i 

(n+l)x2(6-f Nm) 


Xfi and y Fi (i = 1, 2 ...n) represent the coordinates of the i-th node from centre of mass of 
the fuselage. The vector {F} 5 and the matrix [C s ] corresponding to the preselected sensor 
locations are subsets of the vector {Y} and matrix [C] respectively, given above. 
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Definition of condition number 

The quantity j|A|| |[A -1 1| is defined as the Condition number of a nonsingular square 
matrix A (Ref. 41) and will be denoted as Cond (A). In the present study, the condition 
number is defined based on infinity norm of matrix A. 

The infinity norm of a general matrix A of size mx n is defined as 


Mlloo 


max 

Ki<m 
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Appendix C 


Helicopter model as two degrees of freedom system 



The equations of motion of the two degrees of freedom system can be written as 

MgbZgb + K(Zgb ~~ Zp)= Fz + Fc 
MpZp + K (Z F — Zgb)— 0 — Fc 

Writing the equations in matrix form, 


Mob 

0 

\ Zen 


K 

-I< 

f Zgb 

( Fz + F c 


{ .. 

> + 



< 

= { 

0 

M F _ 

{ Zf 

1 

-I< 

I< 

1 J 

l ~ p c J 


Assuming the harmonic input Fz — Fze wt , the response and the control force can be 
written as, 


Z G b = Zcse" 1 , Z F = Z F e M , Fc = F c e<“‘ 
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Substituting the soultion in the equations, the amplitude of the response can be obtained 
from the following relation. 


K - M gb u 2 -K 
-K K- M f cj 2 




F z + F g 1 
~Fc } 


The response can be written as 


Zqb 

1 

K - M f lo 2 

K 


Fz + Fc ] 

J 

MgbMfu 4 — Kuj 2 {Mqb + Mp) 

K 

K - M gb u 2 

< 

i 1 


- (K - M f u) 2 )F z - M f u 2 F c 

GB ~ M gb MM ~ Koj 2 (Mgb + M f ) 


and, 


Thus, 


KF Z + Mgbu 2 Fc 
MqbMfu* — Ku 2 (Mgb Mp) 


Z CB = ZgbF^ 1 
Z F = Zpe™ 1 


The fuselage response Z F will be zero, if 

Pc = ~M^ Pz 

The negative sign in the above equation indicates that F G must be out-of-phase (180 
degrees) to the input force F z 



